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We consider the end state of collapsing null radiation with a string fluid. It is shown that, if 
diffusive transport is assumed for the string, that a naked singularity can form (at least locally). 
The model has the advantage of not being asymptotically flat. We also analyse the case of a 
radiation-string two-fluid and show that a locally naked singularity can result in the collapse of such 
matter. We contrast this model with that of strange quark matter. 
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The end state of a collapsing star has attracted much attention in its own right and especially so since the statement 
of the Cosmic Censorship Conjecture (CCC) (See d0] for a detailed discussion of the CCC). While the end state of a 
sufficiently massive collapsing star is a singularity ]8|, there is no guarantee that such a singularity should be hidden 

■ as is required by the CCC 16] . Indeed, a number of papers have been written describing various scenarios in which 
naked singularities form (See for a recent review.). Such examples clearly violate the CCC. Future formulations 

CN| ■ of the CCC will have to be carefully re-considered in light of these examples. 

Here we present another spherically symmetric spacetime which admits a naked singularity. This spacetime is not 
asymptotically flat and describes the collapse of null radiation with a string fluid first introduced by Glass and Krisch 
. A full physical description of the model is contained in . We merely quote the relevant information for our 
, analysis of the singularities. 

However, we do make one observation: In their analysis Glass and Krisch [(J Q state that the matter content of 
their model is a radiation-string two fluid. The matter content is taken to be a string fluid together with Vaidya 

■ null (pure) radiation. We proceed in a similar manner for the main part of our paper. However, in f3 we consider 
""^ | the special case of Ricci flat spacetimes. Here, we do not assume that the string satisfies a diffusion equation. The 

combination of these two characteristics (Ricci flatness and non-diffusing strings) now allows the matter content to 
be treated as a different radiation-string two fluid (or, perhaps, a three fluid), where, in addition to the Vaidya null 
(pure) radiation, we also have incoherent radiation (as the radiation fluid equation of state p = p/3 is now satisfied) 
together with the string fluid. 

. t-H ! This note is organised as follows: In i[n] we introduce the metric of Glass and Krisch and present the main 
equations to solve. We make a few comments on these equations and provide some special solutions. These solutions 
are used in flllll to study the nature of the collapse. We show that one solution results in a black hole while the other 
could result in a naked singularity. The strength of the naked singularity is examined in illVI During the analysis, we 
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present the general forms of the Kretschmann and Ricci scalars. We investigate the collapse of Ricci flat spacetimes in 
JVl and show that the presence of a string fluid can lead to the occurrence of locally naked singularities for radiation 
fluids. We compare our results to a similar metric that has recently appeared in the literature in WII 

II. THE NULL RADIATION-STRING FLUID MODEL 

Glass and Krisch start with a spherically symmetric metric in retarded time co-ordinate, u, of the form 

ds 2 = Adu 2 + 2dudr - r 2 (di? 2 + sin 2 Mtp 2 ) , (1) 

where 

2m(u,r) , „ 

A = l ^-A (2) 
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Note that they chose the mass function to be a function of both retarded time and distance along the outgoing null 
geodesies. As a result, they were able to extend the Vaidya 0] radiating metric to include a string fluid. 
Remark: A spherically symmetric line element in advanced time u is usually taken to be of the form 

ds 2 = -Adu 2 + 2dudr + r 2 {dd 2 + sin 2 Mip 2 ), (3) 

Thus, the difference between (QJ and is not only a signature change, but also a change of sign of u. We will 
comment further on this later. 

Glass and Krisch Q introduced the unit vectors -0 M ,f M ,i? M and </? M such that 

SV = v^v v - f^fy - d^-d v - tp^ipv, (4) 

where the unit vectors have their usual definition from (T|. The Einstein tensor then took on the form 

2m 2m' „ „ „ m" - - , „ , . 

= ~ l ^ lv t(V» - J > r *) + —yV^Vv + ifi^tpu), (5) 

where ' and ' denote differentiation with respect to r and u respectively and Z M is a principal null geodesic vector, 
from which the nature of the matter content is easily seen. 

Using the Einstein field equations they identified the matter portion of the Einstein tensor as a string fluid: 



Tfj, u = iil^lv + pv^Vy + Prf^fu + p^ifl^v + <Pn0v) (6) 



and thus obtained 



imp = ~ (7) 

m! , , 

Anp = -Anpr = — (8) 

m" 

8lTp± = 9 

r 

where tp is the energy density of the Vaidya radiation, p is the string fluid energy density and p r and p±_ the string 
fluid thermodynamic radial and tangential pressures respectively. 

In addition to the field equations Q - © they also assumed that the density p satisfied the diffusion equation 

P = ^ r [r 2 £), (10) 



r 



where V is the positive coefficient of self-diffusion. Manipulating equations lj7)l-Qll)[l they concluded that, once a 
solution to the linear PDE (fTU|l was obtained, the mass function could be found by integrating © and 



.2 d P 



rh = 4tt2V-^-. (11) 
or 

They presented a few solutions of l|l(J|) and determined the corresponding expressions for m. 

Equation (|llfl is a result of the integrability conditions taking J7J, (jHJ and (|1(J|) into account. However, it would 
seem that an additive function of u was omitted in I jllj l which should be written as (See Appendix |A"1) 

m = AnVr 2 ^ + k(u). (12) 
or 

However, in our calculations, we will still use (|11|) and merely add on an arbitrary fc(tt). 

We observe that l|10f) can be transformed to the heat equation with a nonconstant diffusion coefficient if we set 
R = — l/r. Thus all the known solutions to the heat equation can be applied in this setting. 

It is easy to see that, when p — R(r)U(u) we obtain 
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where C\ and C2 are constants and A is the separation constant. This leads to the following mass function: 



m(u,r) = 

T 



(14) 



If we require that p — R{r) + U(u) then 



with corresponding mass function 



C Q + Xu+^ + ^- (15) 
r 62? 



152? 

If we let 



, \ 27rr 5 A , ^ 2 tt r 2 (2 r(u A + C ) + 3d) „ 
m(u,r) = 42?7r-uCiH — - — — + C 2 . (16) 



M(u) = -AVttCiu (17) 



be the Vaidya mass and note that the 47rCor 3 /3 term makes a contribution similar to the cosmological constant, we 
can define 

s( „, r) = ^ + 2 " 3 ^-" + 3ft> +a (18) 

as a mass contribution from the string fluid. We note that both S(u) and M(u) have a term involving C±. Thus, in 
this model, it is not possible to separate the Vaidya mass from the string mass. This is a result of the relationships 
© and (fTT|) which are a consequence of itTUl) . 

We note that, if we require the string to behave like a perfect fluid, then we must have p r = p± = p and so p = —p 
remains as the equation of state for the string fluid. Using equations © and ©, this means that 

rm" - 2m' = (19) 

whence 

m = F{u)r 3 + G(u). (20) 

If we further impose 1|1(J[) . we have that F(u) — Fq, a constant and the resulting part of the mass function, ie. For 3 is 
a cosmological constant term - the metric is Vaidya-de Sitter (See |2(| for a discussion of the CCC in this context.). 
In the resulting discussion, we will assume that the fluid is anisotropic. 

III. NATURE OF THE COLLAPSE 

The geodesic equations (taking the null condition K a K a into account, where K a — dx a /dk is the tangent vector 
of a geodesic) for the metric are 



= (21) 



dfe r 2 \dfc 

«. 7 (S; =»• < 22 ' 

To investigate the structure of the collapse we need to consider the radial null geodesies defined by ds 2 = taking 
(j) = 6 = into account. For (|TJ the radial null geodesies must satisfy the null condition 

du 2 (23 ) 
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However, we recall that the co-ordinate u is a retarded time co-ordinate. As remarked earlier, in most examples in 
the literature, an advanced time co-ordinate is used. In order to compare our results with those usually found in the 
literature we need to set u — * — u in 1)23(1 and obtain 

*L = I 

r 

For appropriate forms of to, this equation will have a singularity at r — > 0, u — > 0. In order to classify the radial 
and non-radial outgoing non-spacelike geodesies terminating at this singularity in the past, we need to consider the 
limiting values of X — u/r along a singular geodesic as the singularity is approached ^j- Thus, for the geodetic 
tangent to exist uniquely at this point we must have that 

Xq = lim — = lim = lim —^7 7. (25) 

u->0 r^O r ii^O r^O dr u^O r^O I _ 2m(-u,r) >■ ' 

r 

Thus the crucial aspect of metric JQ) that dictates the nature of the collapse is, as expected, the mass function to. 

We note that p = po, a constant, is a solution of ((11) I) . The corresponding mass function is m(u, r) = mo + Air pgr 3 / 3. 
For this form of to, we have 

2 2 

X ° = u 1 2m„+8w 3 73 = J _ Um ^ = °- ( 26 ) 

r r 

Since there is no real, positive Xq, there is no non-spacelike geodesic emanating from the singularity and so the 
singularity is not visible to any observer. However, this is not a realistic model as to is independent of time and so 
the model is static. The reason we comment on this form of m is to show that any constant part of to will result in 
a black hole. As a result in our further analysis of the structure of the collapse we will always set too = 0. 

Let us consider solution 1|14|) . It is clear from that solution that, as r, u — > 0, we have to — > k, a constant and so, 
from ((26(1 that a black hole forms. 

In the case of ijTfil l (with C2 = 0) we observe that l|2"5|) reduces to 



from which we obtain 



Xo = a±= 1±V1 - 64 ^. (28) 
* 16Ci7rX> V ' 



If we define 5 = 8Ci7r2? we can rewrite ((28(1 as 



1± vT" 



Xo = a ± = . (29) 

In order for real positive solutions in l(29|) to exist we must have that 

5 < ~- (30) 

This is exactly the same form of Xq (and restriction) obtained by Dwivedi and Joshi Q and Wagh and Maharaj 
|20|. That the similarity arises should not be of any surprise. It is clear that, if one is dealing with a mass function 
that is of the form Xu + F(u, r) where F(u, r)/r is polynomial in r and/or u (of degree at least one), then the only 
contribution to the condition for the existence of a time-like geodesic emanating from the singularity is from Am, the 
Vaidya mass. Indeed, this is the reason that asymptotic flatness is not essential for the existence of a locally naked 
singularity, a conclusion first obtained in |20| . 



IV. STRENGTH OF THE SINGULARITY 



The main importance of determining the strength of the singularity is due to the fact that the CCC does not need 
to rule out the possibility of the occurrence of weak naked singularities |l4j . This arises as one may continue the 
geometry of a weak naked singularity through the singularity to make it geodesically complete 0,0] (See also 01-)- 
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In order to determine the strength of the singularity we utilise the ideas of Nolan 0] as explained in Thus, for 
a weak singularity, we require 



Using the definition 



we have, asymptotically, 



dr 

— ~ d => r ~ d fc- (31) 



7/ (17/ 

X Q = lim - = lim ^ (32) 



— ~ d X => u ~ d X Q k. (33) 



The geodesic equation <|21ll then becomes 



cgw ^ 2Xg7Td / 4Ad 4 fc4 + 20^fc 3 dX A + 20d 2 r.k 2 VC a + 15d kVd + 30d 2 X a C 1 ) (34) 
dk z lbkV 

which is of 0(k)~ 1 and so is inconsistent with 1)33(1 . Thus the singularity is gravitationally strong in the sense of 
Tipler [13. 

The Kretschmann scalar K = R a bcd,R abcd for (JIJ is given by 

K = -jr((m") 2 ~ 4m"m'r 3 + 4mm" r 2 + 8(m')V - Wmm'r + 12m 2 ). (35) 
If we utilise the mass function (|16|l . I|35(l evaluates to 

K = 64 ^ 2 6 (1200r 6 X> 2 MAC + 600r 6 X>VA 2 + 900r 5 X> 2 C C*i + 300r 8 A 2 X>u + 300r 8 AX>C + 150r 7 AX>Ci 
+ 2700X>VC 2 + 720r 5 AX> 2 uC*i + 600r 6 X> 2 C 2 + 450r 4 X> 2 C 2 + 53r 10 A 2 ) 

which clearly diverges at the naked singularity. Hence the singularity is a scalar polynomial singularity. 
The Ricci scalar for Q is given by 

R=^-(m"r + 2m') (36) 



which, in the case of (|lfc>|) evaluates to 



R = (r 3 A + AVruX + AVrC a + 3X>Ci) (37) 

Vr 



and also diverges at the naked singularity. 



RICCI FLAT SPACETIMES 



It is clear from <|36|) that none of the solutions presented above cover the case of Ricci flat spacetimes. If we relax 
the requirement that the string diffuses (ie. we do not impose (|1U|) ^ we can now consider the case of R = 0. From 
(I3GII this means that the mass function must take on the form 

m(u,r)=M(u) + ^-, (38) 
r 

where we can again interpret M(u) as the Vaidya mass and S(u) as the mass contribution from the matter content. 
Forcing the spacetime to be Ricci flat requires the matter content to obey the (incoherent) radiation fluid equation 
of state, ie. 

_ p r + 2 P± Pr P 

= -y = 3 (39) 
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TABLE I: Positive numerical solutions to I4H for specific forms of a = 1.0 and /3. 



P X 

-1 1.35321 

-0.9 1.42537 

-0.8 1.51454 

-0.7 1.62861 

-0.6 1.78145 

-0.5 2.0 

-0.4 2.3436 

-0.3 2.96715 

-0.2 4.3797 

-0.1 9.1457 



from JSJ and @ . Thus if we begin with a metric of the form Q (together with an energy-momentum tensor of the 
form <©) and require it to be Ricci fiat, then, as expected, the fluid must satisfy (|39|1 ■ These spacetimes must be 
asymptotically flat via (|38|) . 

To proceed further, we must assume some form for M(u) and S(u). Given the interpretation above, it makes sense 
to assume that M(u) = —cm. In order for S(u) to contribute to the existence (or otherwise) of a naked singularity, 
it must be of 0{u) 2 . Thus we set S{u) = /3u 2 . The null condition it%4l now becomes 

du _ 2 

— ^ 2mi+2(3u 2 /r ' \ ' 

r 

This equation has a singularity at r — > 0, u — > 0. The condition for a geodetic tangent to exist uniquely at this point 
reduces to 

u du 2 2 
An = hm — = hm — = hm „ , = ^ (41) 

u^o r^o r u^o r^o dr u^o r~>o i _ 2 ""+ 2 / 3 " 2 A 1 - 2aX - 2(3X?; 

r 

This cubic equation in X has (at least) one real positive root provided that, for a > 0, (3 < 0. This can be seen from 
the numerical results [2lJ given in Table H| Thus, in the case of Ricci flat spacetimes, we also have the possibility of 
a locally naked singularity. 

As above, we investigate the strength of the naked singularity. Taking l|31fl and (|33fl into account, (|21|l (with i|38|) 
and the assumed forms of M(u) and S(u)) becomes 

d^ 2(3X - a 

dk 2 d k [ ' 

which is inconsistent with \'.\'M and so this singularity is also gravitationally strong. 
For Ricci flat spacetimes, the Kretschmann scalar l|35|) is given by 

K = {US(u) 2 + \2S{u)M{u)r + 3M(u) 2 r 2 ) (43) 

which, for the assumed forms of M(u) and S(u), evaluates to 

16u 2 

K = —^{Ufl 2 u 2 - 12f3aur + 3aV). (44) 
Clearly l|44|) diverges at the naked singularity. As expected, we have a scalar polynomial singularity. 



VI. DISCUSSION 



We have presented a scenario for the gravitational collapse of a string fluid together with Vaidya null radiation in 
which the end state becomes a naked singularity. We have also investigated the collapse of Ricci flat spacetimes. Such 
collapse can also give rise to naked singularities. Both scenarios provided naked singularities which were gravitationally 
strong and, as expected, were scalar polynomial singularities. These are two more in the long list of examples that 
proliferate the literature and should serve to tighten future formulations of the CCC. 
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We comment briefly on our results and those of Ghosh and Dadhich 5]. The metrics (apart from a signature 
change and change of sign of time co-ordinate) are similar. As a result, the Einstein field equations also reduce to 
(Q-EJ) (though their p represents isotropic pressure). However, in their model, Ghosh and Dadhich Q interpret the 
matter content as strange quark matter. Thus, while ip is still the radiation density, their p and p refer to strange 
quark matter thermodynamic pressure and energy density respectively. As a result of that interpretation, they impose 
an equation of state relating p to p: 



where B is the bag constant. This differential equation can be treated as an ODE and the general solution was 
presented in 

Here, following Glass and Krisch |(|, the matter content has been interpreted as a string fluid which diffuses like 
point particle diffusion where the number density, n, diffuses from higher numbers to lower according to llOfl (Note 
we have used the fact that p — M s n, M s a constant, to obtain Q10[l.). The string equation of state is always satisfied 
as a result of the energy- momentum tensor (0. As a result, the solutions presented here (and in 0) will not coincide 
with those found in J5|. 

That other interpretations of the metric QJ together with 10 are possible is evident from the case of Ricci flat 
collapse investigated in <J3 Indeed, this case is of considerable importance as the matter content satisfies both the 
string equation of state p r — —p and the (incoherent) radiation fluid equation of state p = p/3 (see 11390). Thus the 
matter content is a two-fluid, in addition to the Vaidya null (pure) radiation. We have shown that a naked singularity 
can arise in this scenario. In 2] it was shown that a naked singularity only forms if the ratio of pressure to density 
was less than —1/3. Here, we have shown that a naked singularity can exist outside this range, ie. for a ratio of 
1/3. We conclude that it is the presence of the string matter that contributes to the occurrence of the locally naked 
singularity. 

As observed in i HTI an arbitrary function of u, k(u) say, can be added to the mass function m. This additive function 
of u can be seen as a further contribution from the string mass S. However, unless it is linear in u, it would not 
contribute to the condition (|3"U)l and even then, would only serve to redefine 5. 

We note that the usual manner in which the existence of a naked singularity is presented is to show that the 
singularity forms before the event horizon forms. Interestingly, it has been shown that, for some models, an apparent 
horizon will not form at all 

0,110 (See also 
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Here we consider the derivation of l|ll|) . Differentiating J7J) and © with respect to r and u respectively results in 



p=-{p-AB), 



(45) 
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APPENDIX A: INTEGRABILITY CONDITIONS FOR m(u,r) 



m, 



m, 



-47r(r 2 V>)' 
4:nr 2 p. 



(Al) 
(A2) 



Thus 



(r 2 V)' 




(A3) 



However, (|1U|) yields 




(A4) 



Thus, from (|A"Tj) . l|A3jl and (fAlJ) we have 




(A5) 
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whence 

(•*»)■ (A6) 

Glass and Krisch [(| presumably integrated this equation to obtain (|llf> . However, since this is a pde, the integration 
should yield 
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